Forced and self-excited oscillations of an optomechanical cavity 
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We experimentally study forced and self-excited oscillations of an optomechanical cavity which is 
formed between a fiber Bragg grating that serves as a static mirror and between a freely suspended 
metallic mechanical resonator that serves as a moving mirror. In the domain of small amplitude 
mechanical oscillations, we find that the optomechanical coupling is manifested as changes in the 
effective resonance frequency, damping rate and cubic nonlinearity of the mechanical resonator. 
Moreover, self-excited oscillations of the micromechanical mirror are observed above a certain optical 
power threshold. A comparison between the experimental results and a theoretical model that we 
have recently derived and analyzed yields a good agreement. The comparison also indicates that the 
dominant optomechanical coupling mechanism is the heating of the metallic mirror due to optical 
absorption. 



I. INTRODUCTION 

Studies combining mechanical elements in optical reso- 
nance cavities Ql experience a significant surge in pop- 
ularity in recent years due to the fast progress made in 
both microelectromechanical systems (MEMS) and opti- 
cal microcavities. For example, optomechanical coupling 
of nanomechanical mirror resonators to optical modes 
of high-finesse cavities mediated by radiation pressure 
has a promise of bringing the mechanical resonators into 
the quantum realm [3-llJ- Furthermore, the micro- 
optoelectromechanical systems (MOEMS) are expected 
to play an increasing role in optical communications fl^l 
and other photonics applications [13-15]. 

In addition to the radiation pressure, another impor- 
tant force that contributes to the optomechanical cou- 
pling in MOEMS is the bolometric force [l6l - [2^ . also 
known as the thermal force. This force can be attributed 
to the thermoelastic deformations of the micromechani- 
cal mirrors. In general, the thermal force plays an im- 
portant role in relatively large mirrors, in which the ther- 
mal relaxation rate is comparable to the mechanical res- 
onance frequency. Phenomena such as mode cooling and 
self-excited oscillations have been shown in systems in 
which this force is dominant [16-1^, [23|. Existing theo- 
retical models that describe these phenomena quantita- 
tively [l^, [l^-Ull are based on energy or harmonic bal- 
ance methods which provide good predictions of the sys- 
tem steady state, but lack the ability to fully describe its 
complex dynamics. 

Recently, we have developed a slow envelope dynamical 
model of an optomechanical system which includes radi- 
ation pressure, thermal force, and changes to mechanical 
frequency due to absorption heating [2g. The theoreti- 
cal predictions, which are derived using a combined har- 
monic balance and averaging method [27], include all the 
experimental phenomena shown by optomechanical sys- 
tems with a bolometric force, such as linear dissipation 
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renormalization and self-excited oscillations. In addition, 
the model enables prediction of additional nonlinear ef- 
fects, namely, the change in the sign of the mechanical 
cubic nonlinear elastic and dissipative terms as function 
of the optical power incident on the cavity and its exact 
detuning from resonance |23J. 

Here, we present experimental results that demon- 
strate all the major dynamical phenomena which are the- 
oretically implied in Ref. [1^. In order to facilitate the 
study of optical cavities with micromechanical mirrors 
spanning a wide range of different geometries and ma- 
terials, we employ a fiber Bragg grating (FBG) [28| as 
a static mirror of the optical cavity. The wave length 
dependent transparency of the FBG allows us to achieve 
different coupling conditions between the optical mode 
inside the cavity and the incident light, thus effectively 
controlling the cavity's finesse. 

A very reasonable fit between theory and experiment 
is achieved using two distinct geometries of the microme- 
chanical mirror composed of two different metals - AuPd 
and aluminum. The fits include changes in the linear dis- 
sipation, the threshold, the frequency and the amplitude 
of self-excited oscillations, and the thermally induced fre- 
quency shifts under different conditions. In addition, we 
show optically induced changes in the nonlinear response 
of the micromechanical mirrors. 



II. EXPERIMENTAL SETUP 

In the investigated system, an optical resonance cavity 
is created between a suspended metallic micromechanical 
mirror, which is free to oscillate in a direction parallel to 
the optical axis, and a stationary mirror in the form of 
a FBG as shown in Fig. [1] The system is located in a 
vacuum chamber inside a cryostat with a typical pressure 
of 3pbar and temperature of 77 K. 

A micromechanical mirror is fabricated on a silicon- 
nitride membrane using electron beam lithography and 
thermal evaporation of metal. Following these steps, the 
membrane is removed by electron cyclotron resonance 
(ECR) plasma etching, and the micromechanical mirror 
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FIG. 1. (Color online) (a) The experimental system. In ad- 
dition to the network analyzer, a number of measuring and 
excitation schemes can be used, such as homodyne vibration 
detection with a lock-in amplifier, and direct time sampling of 
the reflected optical power. A tunable infra-red laser is used, 
with a maximum output power of 2mW. A single mode op- 
tical fiber is used to transmit light to and from the sample, 
allowing reflection measurements (with the aid of a circulator 
to separate the incident and the reflected beams). A graded 
index lens at the end of the optical fiber is used to focus the 
light on the mirror. The focus distance is ~ 40 pm. A FBG 
located inside the fiber can serve as a second mirror to create 
a relatively high-finesse cavity at wavelengths that fall inside 
the Bragg region. The micromechanical mirror can be ca- 
pacitively actuated by applying a voltage between the mirror 
itself and a ground plate located 500 ym below it (the ground 
plate is not shown). Panels (b) and (c) exhibit the top views 
of the micromechanical mirrors employed in the experiments. 
The thickness of the metal layers is 300 nm for AuPd samples, 
and 200 nm for Al samples. 



becomes suspended. This fabrication process is similar 
to the one described in (29| . 

Two main suspended mirror configurations were used 
in our experiments, a gold-palladium (Auo.85Pdo.15) rect- 
angular mirror and an aluminum doubly clamped wide 
beam. The dimensions of the devices are given in Fig. [T] 

The micromechanical mirror can be actuated capac- 
itively in the direction of the optical axis by applying 
voltage between the mirror and the ground plate of the 
package used to mount the sample in the vacuum cham- 
ber. The ground plate is parallel to the mirror and lo- 
cated 500 Jim below it. 



The FBG, which has a length of Lb = 20 mm, is 
formed using a phase mask having a period of 1062 nm 
on a single mode optical fiber having effective refractive 
index of rioff — 1.46 near the wavelength of 1550 nm. A 
microlens made from a section of a graded index fiber 
having length of 0.45 pitch is spliced to the the end of 
the fiber [30]. A brief analysis of the FBG's optical prop- 
erties is given in Sec. IIII Al 

The optical fiber can be moved in three orthogonal di- 
rections by the means of piezomotors with an accuracy 
of approximately 1 nm. Generally, the fiber is positioned 
above the center of the micromechanical mirror at a focal 
distance of the microlens, which is « 40 pm. The length 
of the optical cavity can be changed by moving the fiber 
along the optical axis. We control the wavelength and 
the power of the light incident on the cavity by using a 
variable wavelength infra-red laser and a variable fiber- 
optic attenuator, respectively. The light reflected off the 
cavity back into the fiber is separated by the means of a 
circulator and converted to an electrical signal by a pho- 
todetector. The experimental system is shown in Fig. [TJ 

The finesse of the optical cavities created in the pre- 
sented experiments is of order of ten. We estimate the 
optical relaxation time to be of order 10~^^ sec. It fol- 
lows that optical retardation can be neglected in our sys- 
tem, and thus the optical energy stored in the cavity is a 
function of the momentary displacement of the microme- 
chanical mirror. 



III. THEORETICAL MODEL 

An extensive theoretical analysis of the dynamics of 
a micromechanical oscillator acting as a mirror in a low 
finesse optical cavity based on a slow envelope approxi- 
mation can be found in Ref. [2^ . Here, we state the main 
results from that work, and present a short discussion on 
the FBG optical properties. 

A. Optical cavity 

The finesse of the optical cavity is limited by loss mech- 
anisms that give rise to optical energy leaking out of the 
cavity. The main escape routes are through the FBG, 
through absorption by the metallic mirror, and through 
radiation, and the corresponding transmission probabil- 
ities are respectively denoted by Tb, Ta and Tr. The 
transmission probability Tg through the FBG is evalu- 
ated using the coupled mode theory (28l . l3l| 



where ds = {oj — ujb) LbtIcs/c is the normalized detun- 
ing factor, oj and ojb are respectively the laser and Bragg 
angular frequencies, c is light velocity in vacuum, and Vb 
is the FBG coupling constant. 
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Let a; — a:o be the displacement of the mirror relative 
to the point xq, at which the energy stored in the optical 
cavity in steady state obtains a local maximum. For a 
fixed X the cavity reflection probability Rc, i.e. the ratio 
between the reflected (outgoing) and injected (incoming) 
optical powers in the fiber, is given by 



In our experiments, however, the contributions of purely 
mechanical nonlinearities are negligible, as will be shown 
in Sec. IVC] 

Following Ref . [1^ , we write the equation of motion as 
x+'^x+uj'^x = 2fjnCOs{ujo+ao)t+F,.p{x)+Fth{x), (5) 



Rc 



Tb-Ta-Tr 
Tb+Ta+Tr 



1 — cos 27r- 



1 + 2 (^) (l-cos27r^; 



(2) 



where L is the distance between two successive resonance 
positions of the micromechanical mirror (i.e., half the 
wavelength), and 



T^{Tb+Ta + Tr) 



L 
2^ 



is the full width at half maximum parameter. The ef- 
fective optical power I{x) impinging on the suspended 
micromechanical mirror can be expressed as 



I{x) = 



Imax ( 2 ) 



where 



cos 27r2^] 



Cfc^pump; 



(3) 



is the maximum optical power incident on the mirror, 
-^pump is the power of the monochromatic laser light in- 
cident on the cavity, and 



{Tb + Ta + TrY 



is the resonant enhancement factor of the intra-cavity 
power. Note that for the case of critical coupling, i.e., 
the case where Tb =Ta+Tr, C^e = L/ttT. 

The optical power I(x) is a periodic function, which 
can be approximated by a truncated Fourier series 



I{x) 



E 



Cke 



(4) 



where ftmax should be of order of the finesse or larger 
for the truncation error to be negligible. As shown 
in Ref. [ H, Ck = / maxX«''''e- J^'^''="°/-^, wh ere a = 
1 + h- v'(l + /i)2-l, X = h/^{l + h)^-l, and h = 
^2r2/2L2. 



where a dot denotes differentiation with respect to time 
t, ujQ is the mechanical resonance frequency of the mir- 
ror, Q is the mechanical quality factor, ujm is the tem- 
perature dependent momentary resonance frequency, 
is the external excitation force, and ctq is a small de- 
tuning of the external excitation frequency from wq, i.e., 
ctq "C Wo- The forces resulting from coupling to an opti- 
cal resonance cavity are the radiation pressure Fy-p, and 
Fth, which is a thermal force that appears due to tem- 
perature dependent deformation of the micromechanical 
mirror [si - fSTj . 

In a wide range of micromechanical resonators, inter- 
nal tension can strongly affect the resonance frequencies 
[38l.[39|. Such systems include the doubly clamped beams 
and rectangular mirrors with four suspension beams used 
in our experiments. Changes in the temperature of such 
devices result in thermal expansion or contraction, which 
in turn cause changes in internal tension. These changes 
give rise to a strong temperature dependence of the me- 
chanical resonance frequencies, as will be shown in SecfVl 
For small temperature changes, the momentary mechan- 
ical frequency w„j is assumed to be linearly dependent on 
the temperature: 



= Wo - /3{T - To), 



(6) 



where /3 is a proportionality coefficient, T is the effective 
temperature of the mechanical oscillator, and To is the 
temperature of the supporting substrate. In our samples, 
a significant pretension exists due to thermal evaporation 
process used to deposit the metals during the manufac- 
turing [33', ^sS, 'io'l . The pretension is further increased by 
cooling the samples to 77 K. It follows, therefore, that (3 
is positive in our experiments, i.e., heating of the sample 
reduces its resonance frequency. 

The effective temperature changes can be described by 
the following equation. 



^K{T~To)+vIix), 



(7) 



where k is the effective thermal conductance, and 77 is 
the effective radiation absorption coefficient. The formal 
solution of Eq. (O can be shown to be 



B. Equations of motion 

The micromechanical mirror can be approximately de- 
scribed as a harmonic oscillator with a single degree of 
freedom x operating near primary resonance, which is 
subject to several forces arising from coupling to the op- 
tical resonance cavity. In general, a stand alone microme- 
chanical resonator can exhibit nonlinear behavior l32ll33ll. 



T-To = 77 / /(x)e«(^-*)dr. 



where the initial transient response term 
g-Kt jjij^^ = 0) — To] has been dropped as insignifi- 
cant to the long timescale dynamics of the system. This 
integral relation can be further simplified using the 
slow envelope approximation. The reader is referred to 
Ref. |2a| for further details. 
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Finally, we introduce the radiation dependent forces. 
The radiation pressure force can be expressed as 

F,pix) = vl{x), (8) 

where 

2 

mc 

and where m is the effective mass of the micromechan- 
ical mirror. Light absorption by the mirror has been 
neglected. The thermal force Fth is assumed to be linear 
in the temperature difference T — Tq, i.e., 

Fth = 0{T - To), (9) 

where is a coefficient of proportionality. 

The numerical values of all the physical constants in- 
troduced above will be evaluated in Sec. IIVI 



where cj) = <j) — a^t (the detuning ctq is assumed small) , 

(37] 

n = (jjo Po=wo~Awo, (12) 

K 

and 

where Jn{z) is the Bessel function of order n. 

The term Aujq represents a small mechanical frequency 
correction due to the averaged heating of the microme- 
chanical mirror vibrating with an amplitude Ai. As will 
be shown in Sec.|Vl this correction accounts for the dom- 
inant part of the resonance frequency shift measured in 
our experiments. 



C. Slow envelope approximation 



D. Small amplitude oscillations 



Following Ref. [26| , the dynamics of the micromechan- 
ical mirror can be approximated by a harmonic motion 
with slow varying amplitude and phase, i.e.. 



where 



x{t) ~ Aq + Ai cos ip, 



(10a) 



(10b) 



and where Ai and (f> s^re the oscillator's amplitude and 
phase dll, respectively, and Aq is the static mirror dis- 
placement due to the action of the radiation dependent 
forces. 

By introducing Eqs. (|10p into Eq. ([5]) and using a com- 
bined harmonic balance - averaging method (27| , one can 
derive the following relations which describe the slow 
envelope behavior of the mirror [see Eqs. (22)-(25) in 
Ref. f26l|]: 



Ao 



Ai 



1 



CiJqK 



jAi +Po\v 



(11a) 



Pi-n 



LOq 



k2 -I- Ujf 



2(3 Ao + — 

\ ^0 



■ sm ( 



Wo 



(lib) 



and 



(To + Awo + P2/3?7 



2 

K'^ + UJq 



2(3 Ao ■ 



ojo 



A, 



ujo 

— COS0, (11c) 
Wo 



The evolution equations (jlip can be conveniently sim- 
plified if the vibration amplitude of the micromechanical 
mirror is small compared to the optical resonance width 
parameter F. In this case, 



lo f Or] 



and 



A 



Is 



ojo 



(13a) 



1 .3 /" 



Ais((} = - (ao + Aujs) Au + -^Ai, - ^ cos (13b) 

4 Wo 

where the subscript 's' denotes small amplitude oscilla- 
tions, lo — I{x — 0), and where 



Awos = — /o, 



(14a) 
(14b) 



7=^+^/^^ /^^o. + -^)/^, (14c) 
2Q -I- Wq \ 2wo/ 



AWs = Awqs 



and 



T- + 2 2 /^^o. + TT- 

ZUJo -f Wq \ 2wo 



(3r]3K^ + &4 
q — —- : — —ii 



2k k2 + 4w2 



(14d) 

(15a) 
(15b) 
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Note that a prime denotes differentiation with respect to 
X, i.e., Iq = dl(x = 0)/dx. 

The evolution equations ([T^ describe a Duffing-hke 
nonhnear oscillator with nonlinear damping (ssl . l4l| In- 
terestingly enough, the sign of the nonlinearities depends 
on the sign of the second derivative of Iq with respect 
to X. For example, Eq. (|15ap predicts that the system 
should exhibit hardening behavior near the maximum of 
the optical resonance (more precisely, in the region where 
Iq < 0, i.e., |a;o| < T/2^/3), and softening behavior oth- 
erwise. This effect is experimentally illustrated in Sec. |V] 
for both types of micromechanical mirrors studied. 

Another interesting effect that depends on the optical 
detuning xq of the micromechanical mirror is the change 
in the effective linear damping coefficient 7 as function 
of Iq, which is evident from Eq. (|14c|) . From the exper- 
imental point of view, it is convenient to introduce the 
effective quality factor as 



1 



icS 



27 



(16) 



We expect an increase in l/Qoff as compared to the 
purely mechanical value 1/Q in the region in which /g > 
(xq > 0). corresponding to the mode "cooling" effect 
d, H, [21I, \^ , and, conversely, decrease in the effective 
dissipation for the values of xq at which Iq < 0, i.e., xq < 
0. In this region, the effective dissipation may become 
arbitrarily small and even change sign, resulting in a Hopf 
bifurcation followed by possible self-excited limit cycle 
oscillations |i26i. i43i] . 



E. Self-excited oscillations 

Self-excited oscillations may occur in a swtem de- 
scribed by Eqs. (|lip if a stable limit cycle \I§L ex- 
ists in absence of external excitation. In other words, a 
nonzero solution of the following equation, together with 
Eq. piap . is required: 



2Q ^ k2 



Uq 



PiV- 



Uq 



2/3 Ao + — 

UJQ 



= 0. (17) 



Again, we refer the reader to Ref. [26] for a full analy- 
sis of different bifurcations and limit cycles types which 
may appear in systems under study. In our experiments, 
a single stable limit cycle is observed, appearing beyond 
the threshold of a supercritical Hopf bifurcation. As 
expected, the region in which the system develops self- 
excited oscillations coincides with the region of negative 
effective linear dissipation, i.e., 7 < [see Eq. (|14cp ]. in 
which the zero amplitude solution = is unstable. 
For a given nonzero oscillation amplitude Ai, the os- 



cillation frequency correction is given by Eq. ([TTc]): 



I\ 



V^TT^ (2/3^0 + — 



UJQ 



-Aujt. (18) 



IV. PARAMETER EVALUATION 

We now turn to evaluate the physical parameters, 
which are defined in the previous section, for a rectan- 
gular Auo ssPdo 15 mirror, whose dimensions are given in 

Fig.m 

The environment temperature in our experiments is 
To = 77 K. Using a weighted averaging of the values for 
gold [isl . l46l | and palladium [47] , we estimate the values 
of the density p = 18.2 x lO'^kgm"'^, the mass-specific 
heat capacity Cm = 103Jkg^^K~^, and the thermal 
conductivity k = 281Wm^^K^^ for Auo.85Pdo.15 mir- 
ror at 77 K. Although by no means precise, this simple 
averaging method provides a reasonable accuracy in our 
case. 

We take the effective mass of the micromechanical res- 
onator to be the mass of the mirror (the mass of the 
suspension beams is neglected). Using the mirror's di- 
mensions and the density value derived above, we find 
that TO « 20 ng. 

The effective thermal relaxation rate k can be evalu- 
ated as follows: 



300 nm x 5 jim k 



212 pm 



= 3.9 X lO-' 



1 



sec 



In order to estimate the value of the effective radiation 
absorption coefficient 77, the reflectivity of the microme- 
chanical mirror must be known. In the literature ^48i . i49| , 
experimental values between 98% and 99% are given. We 
find that an empirical value of 98.4% fits our experimen- 
tal results. It follows that 



1 - 0.984 



V 



= 7.9 X 10'^ 



K 



The high refiectivity of the micromechanical mirror al- 
lows us to neglect any absorption when estimating the 
radiation pressure coefficient, resulting in: 

2 N 

TOc kg W 

The estimation of the thermal frequency shift coef- 
ficient /3 is not straightforward. The order of magni- 
tude can be estimated by measuring the mechanical res- 
onance frequency of the mirror at room temperature and 
at 77 K (106kHz and 160kHz, respectively), resulting in 
/3 sa 0.012ajo/l K. However, in order to give an accurate 
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estimation of /3 for small temperature changes around 
77 K (or any other ambient temperature), one would re- 
quire preexisting knowledge of the tension inside the sam- 
ple, the exact relation between the tension and the me- 
chanical resonance frequency, and, most importantly, the 
exact temperature distribution inside the sample due to 
nonuniform heating by a focused laser beam. This data is 
not readily available from our measurements. Therefore, 
we treat p as one of the fitting parameters. The best fit 
is achieved for 



O.Olwo 
IK 



1.006 X 10' 



rad 
sec K ' 



which is remarkably similar to the value estimated above. 

Although the majority of the parameters defined in 
Sec. mil can be evaluated using general physical consid- 
erations or direct measurements, 9 is not easily deter- 
mined, because the physical processes responsible for the 
appearance of the thermal force Fth are not well identi- 
fied. Therefore, we derive the value of 6 from experiment. 
The best fit is achieved when 9 is taken to be 



9 = 740' 
By estimating the ratio 



N 



kgK 



10^ 



it follows from Eqs. (fTT|) that the radiation pressure ef- 
fects in our system are negligible compared to the effects 
of the thermal force Fth • 



RESULTS AND DISCUSSION 



sample 


parameter 


value 


Sample I 




2n X 160.088 kHz 


AuPd mirror 


Q 


2.43 X 10^ 


Sample II 




27r X 148.495 kHz 


AuPd mirror 


Q 


6.75 X 10^ 


Sample III 




2-K X 61.25 kHz 


Al beam 


Q 


1 X 10^ 



TABLE I. Mechanical properties of the samples used in this 
study. 



70 nm < r < 140 nm. In general, each time the cavity 
is optically tuned by realignment of the fiber, a slightly 
different finesse can be expected, due to inaccuracies in 
the fiber positioning. 

Instantaneous changes in the micromechanical beam 
displacement x cause changes in the reflected power ac- 
cording to Eq. The signal at the output of the pho- 
todetector can be translated into actual displacement val- 
ues using the calibration discussed above. An example of 
the reflected optical power vs. the optical cavity detuning 
xq, together with sample time traces of mirror oscillatory 
movement is shown in Fig.[2j It is evident from this figure 
that the theory presented in Sec. IIII Al provides a good 
analytical description of the experimental measurements 
of the optical cavity behavior. 



Here, we present a comparison between the experimen- 
tal behavior of three different micromechanical mirrors 
and the theoretical predictions given in Sec. lIIII For con- 
venience, the main mechanical properties of these mirrors 
are summarized in Table ID 

A. Optical resonance cavity 

In our experiments, we tune the optical wavelength 
to a value at which the refiection from the cavity be- 
comes virtually zero at the resonance, a condition known 
as critical coupling. In general, this critical coupling 
wavelength is at the edge of the Bragg region, where 
the FBG reflectivity changes from almost zero to almost 
unity. For example, the optical wavelength used for mea- 
surements of square AuPd micromechanical mirrors is 
1548.83 nm. It follows that the distance between the sub- 
sequent minimums in the reflection Rc (or, conversely, 
peaks in I{x)) is L = 774.4 nm, allowing us to calibrate 
the vertical displacement of the fiber at any temperature. 
A typical finesse of the cavity is between 6 and 11, i.e.. 



B. Linear damping 

We begin our experimental study with investigation of 
what is arguably the most important prediction of the 
theoretical model - the possibility of a significant change 
in the effective dissipation in the vicinity of an optical res- 
onance. In order to measure the effective quality factor 
Qeff defined in Eq. ()16|) at different optical powers /pump 
and cavity detunings xq, we capacitively excite the mi- 
cromechanical mirror at its apparent resonance frequency 
for a short period of time, and then allow the system to 
decay freely to the zero amplitude steady state. During 
this free ring down process, the slow envelope of the me- 
chanical oscillations is measured by the means of a lock-in 
amplifier. The resulting slow envelope is fitted to an ex- 
ponential decay function proportional to e"^'''*, providing 
an estimate of the linear dissipation constant. It is im- 
portant to keep the vibration amplitude small compared 
to r, so the nonlinearities introduced by the detection 
system and the optomechanical coupling [see Eqs. ([TSl) ] 
remain negligible. 

The results presented in Fig. [3] show a good match be- 
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FIG. 2. (Color online) Optical reflection in the vicinity of 
an optical resonance. The incident optical power is /pump = 
3 /iW. Blue circles denote the measured reflected power. Solid 
black line is a theoretical reflection fit [see Eq. ((2)1]. Dashed 
red line represents the theoretical values of the optical power 
/o incident on the micromechanical mirror [see Eq. (13])[. In 
this case, T = 134.5 nm, and the finesse is 5.8. In the inset, 
two time domain traces of the photodetector output signal 
are presented, corresponding to a steady oscillation of the 
micromechanical mirror with a large (thin black line) and a 
small (bold blue line) amplitude. 



tween the experimental values of Qeff and the theoretical 
predictions. The measurement was done at 77 K using 
Sample I (see Table The values of all the system 
parameters used in the fit are similar to those given in 
SecHVl 



C. Nonlinear stiffness and damping effects 




FIG. 3. (Color online) Changes in the effective quality fac- 
tor Qoff [see Eq. (|16ll [ as function of the optical detuning a;o 



and the optical power incident on the cavity, Jp 



Panel 



(a) shows the theoretical value of Q/Qeff, whereas panel (b) 
shows the measured data. In the case of a rectangular AuPd 
mirror presented here, cjo = 27r x 160.088 kHz, Q = 2.43 x 10^, 
and the cavity finesse is 7.3. Panels (c), (d) and (e) show 
cross sections of the top color maps at different values of 
/pump. Blue dots represent the experimental values of Q/Qeff, 
while solid black lines represent the theoretical results. The 
values of all system parameters used in the fit are similar to 
those given in Sec. IIVI The temperature is 77 K. 



The nonlinear effects described in Eqs. psp have been 
observed in all our samples. Here, we present the small 
amplitude frequency response of Samples II and III (see 
Table m, taken at different cavity detuning values. It 
follows from the discussion in Sec. IIIIDI that the elas- 
tic nonlinearity coefficient q should change sign at xq — 
±r/2V3 = ±0.289r. Outside this region the system is 
expected to behave as a softening Duffing-like oscillator, 
while in the region around the optical resonance the be- 
havior should be hardening. The experimental results 
presented in Fig. 2] confirm this prediction qualitatively. 

In general, nonlinear elastic and dissipative effects in 
micromechanical systems can have a non negligible im- 
pact on the dynamics of these systems 132, 33, 41, 50, 51|. 
In the theoretical treatment in Ref. [26|, the nonlinear 
effects which do not stem from optomechanical coupling 
are described by the cubic nonlinearity coefficients aa 
and 73 [see Eq. (7) in Ref. [1^1 . However, the experi- 
mental results show that in our samples the nonlineari- 



ties introduced by the optomechanical coupling are much 
stronger than any preexisting nonlinear effects, at least at 
relatively high optical powers. Therefore, in the present 
work, we have neglected all nonlinearities that do not 
arise due to the interaction with the optical system, i.e., 
as = and 73 — 0. 

D. Self-excited oscillations 

All the samples used in our experiments exhibit the 
phenomenon of self sustained oscillations (i.e., stable 
limit cycle) above certain threshold of the incident op- 
tical power. As expected, these self oscillations always 
occur when < 0, i.e., in the region in which Jg < 0. 
The onset of the self oscillation can be predicted by cal- 
culating the effective linear dissipation coefficient 7 given 
in Eq. p4cp . Self oscillations occur when 7 becomes neg- 
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FIG. 4. (Color online) Changes in the small amplitude non- 
linear behavior as a function of the optical detuning a;o. 
It follows from Eqs. (|15|l that the micromechanical mirror 
should exhibit hardening (q > 0) behavior [Zlj in the re- 
gion |a;o| < r/2\/3 — 0.289r, and softening behavior outside 
this region. This effect is illustrated for two different sam- 
ples, (a), (b) AuPd rectangular mirror measurements, (c), 
(d) Al doubly clamped beam mirror measurements. In both 
cases, the micromechanical mirrors are excited capacitively 
and the frequency response is measured when sweeping the 
excitation frequency up (solid blue line) and down (dashed 
black line). As expected, the frequency responses (a) and (c), 
for which xg > r/2\/3, exhibit softening elastic nonlinearity, 
while the frequency responses (b) and (d), corresponding to 
Xo < r/2\/3, show hardening. The optical power incident on 
the cavity is UfiW for AuPd rectangular mirror and 2.5 ^tW 
for Al doubly clamped beam mirror. 
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FIG. 5. (Color online) The self oscillation amplitude Alc 
as function of the optical detuning xo and optical power in- 
cident on the cavity, /pump. Panel (a) shows theoretical fit, 
whereas panel (b) shows measured data. The color represents 
the value Alc/F. The dashed bold green line in the experi- 
mental color map (b) represents the theoretical threshold of 
self-oscillations. In the case of a rectangular AuPd mirror 
presented here, ujo = 2n x 160.088 kHz, Q = 2.43 x 10^ and 
the cavity finesse is 5.9. Panels (c), (d) and (e) show cross 
sections of the top color maps at different values of /pump. 
Blue dots represent the experimental values of ^lc/F, while 
solid black lines represent the theoretical results. The values 
of all the system parameters used in the fit are similar to those 
given in Sec. IIVI The temperature is 77 K. 



ative. The amplitude and the frequency of the stable 
limit cycle can be found by solving Eqs. ([TT)) and ([T8| . 
respectively. A comparison between the experimentally 
measured self oscillation amplitudes of Sample I (see Ta- 
ble HI and the corresponding solutions of Eq. ([T7|) is 
shown in Fig. [SJ 

It should be emphasized that the theoretical predic- 
tions presented in Figs. [3] and [5] are both based on the 
same set of physical parameters presented in Sec. IIVI and 
on the mechanical properties of Sample I, given in Ta- 
ble U and differ only in the value of F, which changes 
between different experiments, as explained above. It 
follows, therefore, that the theoretical model presented 
here can successfully describe both small vibration be- 
havior and self oscillations with large amplitudes. The 
parameters extracted from experiments in one of these 
two modes of operation can be used to predict the dy- 



namics of the system in the other mode. 

While the theoretical fit shown in Fig. [S] is very rea- 
sonable, a hysteresis phenomenon exists in the experi- 
mental system which can not be explained by the model 
described above. The data presented in Fig.[S]was taken 
while sweeping the optical power from low to high val- 
ues for a fixed value of xq. However, when the optical 
power is swept in the opposite direction, i.e., from high 
to low, the self oscillations disappear at lower values of 
/max. The difference in the threshold optical power can 
be as large as 30%. It should be mentioned that a the- 
oretical analysis of this specific system with parameters 
derived in Sec. IIVI does not predict other stable limit 
cycles, although multiple stable limit cycles (3. fisl [23|. 
as well as subcritical Hopf bifurcations [1^ are possible 
in systems of this type. In the system considered, the 
hysteresis can be possibly attributed to changes in the 



heating pattern and the temperature distribution in the 
vibrating mirror, which cannot be captured by a model 
with a single degree of freedom used in our analysis. A 
multi-mode continuum mechanics analysis of the inves- 
tigated optomechanical system may provide additional 
insight into this hysteresis phenomenon. 

It remains to determine whether the theoretical fre- 
quency shift correctly predicts the corresponding exper- 
imental results both in the case of small vibrations [see 
Eq. (jl4d|) ] and in the case of self-excited oscillations [see 
Eq. HH])]. To this end, we employ Sample II (see Ta- 
ble H]) and measure the spectral power density of the 
reflected light at the vicinity of the sample's mechani- 
cal resonance frequency 148.495 kHz. In this experiment, 
the sample is not excited externally. The incident optical 
power is tuned so the system is expected to develop self 
oscillations at some region of negative optical detunings 
xo- For other values of xo, thermal vibrations manifest 
themselves as a thermal peak, whose frequency is shifted 
by — Aws from ujq. By taking the spectrum traces at 
different values of xq, we are able to measure both the 
frequency of the small oscillations (i.e., the frequency of 
the thermal peak) and the self oscillation frequency. 

The experimental results together with a theoretical 
frequency shift fit are presented in Fig. [S] A very rea- 
sonable fit between theory and experiment is seen in this 
figure. Interestingly enough, the thermal frequency cor- 
rection AojQ [see Eq. ([T^ ] constitutes at least 98% of the 
frequency shift in the entire measured region. 



VI. SUMMARY 

In this work, we experimentally investigate the dynam- 
ics of a metallic micromechanical mirror which is one of 
the two mirrors that form an optical resonance cavity. 
The other, static mirror is implemented as a fiber Bragg 
grating. This unique design allows one to tune the op- 
tical cavity operating conditions to the critical coupling 
domain simply by controlling the wavelength of the inci- 
dent light. 

The finesse of our experimental optical cavities is of 
order ten. Therefore, all optical retardation effects can 
be neglected, and only thermal retardation can play a 
significant role in the dynamics of the micromechanical 
mirror. A theoretical model describing such a system was 
developed in Ref. |26|]. Here, the main results are stated, 
both for small amplitude forced oscillations and for self 
sustained oscillations. 

Theory predicts that coupling of the micromechanical 
oscillator to an optical cavity will result in changes in its 
effective linear dissipation, nonlinear elastic and dissipa- 
tion constants, and the mechanical resonance frequency. 
Stable limit cycles (i.e., self sustained oscillations) will 
occur if the effective linear dissipation becomes negative. 
In addition, multiple limit cycles may be present under 
certain conditions. Two main optomechanical coupling 
mechanisms are postulated, both intermediated by heat- 
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Spectral power density (logarithmic a.u.) 




xo/r 

FIG. 6. (Color online) Spectral power density of the reflected 
optical power as function of frequency and optical detuning 
Xo- The mechanical resonance frequency is 148.495 kHz. The 
incident optical power is /pump = 140 fiW. The color rep- 
resents the spectral power density in arbitrary logarithmic 
units. The region in which self oscillations occur is denoted 
by dashed red lines. The thermal motion peak can be read- 
ily recognized outside this region. The abundance of addi- 
tional peaks in the self-excited oscillations domain may be 
attributed to the nonlinearity of the detection system and 
to mixing between higher mechanical modes. In the inset, 
the experimental values of frequency shift (blue circles) both 
for thermal peak and for self oscillations are plotted together 
with the theoretical predictions (solid black line) given by 
Eqs. (fTid) and (fTS)) . 

ing. The first is mechanical frequency change due to heat- 
ing, the other is a direct force which is a function of the 
temperature difference between the mirror and the envi- 
ronment (thermal force). The radiation pressure force is 
shown to be negligible in our experiments. 

In the present work all the theoretical predictions men- 
tioned above are validated by the means of micromechan- 
ical mirrors with two very different geometries (rectangu- 
lar mirror with four orthogonal suspensions and a wide 
doubly clamped beam). The majority of the physical pa- 
rameters are derived either from general considerations 
or independent measurements. A very reasonable quanti- 
tative agreement between the linear dissipation changes, 
the self oscillation amplitudes, and the frequency shifts 
are achieved. In addition, the theoretically predicted 
changes in nonlinear behavior are demonstrated for both 
mirror configurations. 

Despite the general success of the theoretical fits of 
the experimental data, it is evident that a simple single 
degree of freedom model cannot explain some of the ob- 
served phenomena, most importantly the exact process 
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that gives rise to the thermal force. Another unexplained 
phenomenon is the optical power threshold hysteresis oc- 
curring in the self oscillation measurements. Both effects 
can be possibly attributed to localized changes in heating 
and temperature distribution, and continuum mechanics 
approach is required in order to model them correctly. 
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